clc;
syms x y z;
v=[(x^2*sin(y)),...
    (y^2)*sin(x*z),...
    x*y*sin(cos(z))];
d=divergence(v,[x,y,z]),
c=curl(v,[x,y,z])